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ABSTRACT 

We show how an initially linear spanwise disturbance of the free-stream 
velocity field is amplified by leading edge bluntness effects and ultimately 
leads to a small amplitude but nonlinear spanwise motion far downstream from 
the edge. This spanwise motion is imposed on the boundary layer flow and 
ultimately causes an order-one change in its profile shape. The modified 
profiles are highly unstable and can support Tollmein-Schlichting wave growth 
well upstream of the theoretical lower branch of the neutral stability curve 
for a Blasius boundary layer. 


1. INTRODUCTION 

Efforts to understand boundary layer transition have been underway for a 
number of years now. The relevant experimental measurements are usually made 
on the boundary layer flows over relatively thin flat plates embedded in nomi- 
nally uniform free streams. The final turbulent state arises from the rela- 
tively slow streamwise growth of initially linear instability waves that gradu- 
ally evolve from the small amplitude unsteadiness that either occurs naturally 
in the experiments or can be artificially imposed by the experimenter in order 
to organize the motion. 

An important unresolved issue involves the observed instability wave 
growth that occurs upstream of the theoretical (i.e., two-dimensional flat 
plate) lower branch of the neutral stability curve in some experiments. One 
purpose of the present paper is to provide a possible explanation of this phe- 
nomena. We show that small but steady variations in the upstream velocity 


field can produce somewhat larger streamwise vorticity fields within the bound- 
ary layer which can, in turn, produce significant (i.e,, order-one) variations 
in the streamwise boundary layer profiles. These profiles turn out to be non- 
inf lectional and therefore incapable of supporting rapidly growing inviscid 
instabilities with growth rates that scale with the inverse boundary layer 
thickness. However, the local wall shear becomes small at certain spanwise 
locations which tends to move the lower branch of the Tollmein-Schlichting wave 
neutral stability curve upstream of its theoretical flat plate location. The 
streamwise vorticity field can also amplify the Tollmein-Schlichting waves 
through the resonant mechanism suggested by Nayfeh (1981). The resulting 
three-dimensional disturbances would then grow on the inverse spanwise length 
scale of the disturbance which could (numerically if not asymptotically) be 
equal to the Tollmein-Schlichting wave growth rates. 

Goldstein (1983) and Goldstein et al. (1983) showed that long wavelength 
free-stream disturbances can interact with the rapidly developing leading edge 
boundary layer to generate spatially decaying asymptotic eigenmodes that even- 
tually turn into spatially amplified Tollmein-Schlichting waves once they reach 
the lower branch of the neutral stability curve. The present analysis shows 
that small nonuniformities in the oncoming stream can cause the relevant branch 
to move well upstream of its theoretical Blasius location. 

Crow (1966) also considered the steady distortion of flat-plate boundary 
layer flows by free stream nonuniformities, but that work is quite different 
from ours. One important difference is that the present work is concerned with 
the augmentation of instability wave growth while Crow is concerned with the 
development of three dimensionality. Our emphasis is therefore on the linear 
and nonlinear amplification mechanisms that ultimately produce order-one 
changes in the boundary layer flow, while the boundary layer flow is linearized 
about the Blasius flow in the Crow analysis. In this regard, it is worth not- 
ing that leading-edge bluntness effects play a central role in the present 
work, while Crow considers only an infinitely thin flat plate. Another impor- 
tant difference is that the present work qualitatively predicts the destabili- 
zation of the boundary layer flow, while Crow's analysis, which is only con- 
cerned with the three dimensionality, does not. 

We assume that the dimensions of the 'leading edge ellipse' are of the 
order of the spanwise length scale, say X, of the upstream disturbance field 
and that the Reynolds number based on X, say Rx, is large. Then the upstream 
distortion interacts linearly with the leading edge with the resulting flow 
being well described by the usual 'rapid distortion' theory (Hunt and 
Carruthers, 1990, and Goldstein, 1978). The relevant analysis was, at least in 
principle, given by Lighthill (1956) who showed that the upstream distortion 
produces a cross flow velocity field that becomes logarithmically infinite at 
the surface of the plate. The logarithmic singularity must ultimately be 
removed by viscous effects, which (as in Toomre, 1960) are confined to the vis- 
cous boundary layer region (with the Reynolds number - amplitude scaling being 
considered herein). Our analysis shows that the inviscid cross flow effects 
produce only a linear perturbation of the boundary layer flow in the vicinity 
of the leading edge, where the undisturbed boundary layer undergoes its most 
rapid streamwise development, but that they produce an order-one change in the 
mean boundary layer profiles at large distances downstream where its streamwise 
development is on a considerably slower scale. 
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However, the linear rapid distortion theory solution, which provides an 
adequate description of the external inviscid flow in the vicinity of the lead- 
ing edge, begins to break down at large streamwise distances with the breakdown 
moving further upstream as the surface of the plate is approached. A new non- 
linear solution then has to be obtained in order to describe the inviscid flow 
outside the boundary layer in the physically interesting region where cross 
flow effects produce a significant (i.e., order-one) profile change. The 
thickness of this nonlinear inviscid region is small compared to its streamwise 
dimension but large compared to the boundary layer thickness. It serves as a 
kind of 'blending layer' that connects the boundary layer solution to the lin- 
ear rapid distortion theory solution, which applies at an order-one (on the 
scale of X) distance from the wall. 

The blending layer flow is governed by the nonlinear 'water wave' equation 
whose solution eventually develops a singularity at a certain spanwise location 
and at a finite downstream position due to the well-known wave steeping effects 
associated with that solution. This also produces a singularity (signified by 
the vanishing of the wall shear) in the boundary layer flow at the same (span- 
wise and streamwise) location. The boundary layer therefore, develops local- 
ized regions of relatively small spanwise extent in which the flow is destabi- 
lized as a result of large reductions in the wall shear. 

New local solutions to the blending and boundary layer problem have to be 
worked out in order to understand the structure of these regions. The two 
solutions must then be matched in an appropriate overlap domain. The boundary 
layer solution actually develops a double layer structure as it approaches the 
singularity so that the overall asymptotic structure has three layers in the 
vicinity of the singularity. A second purpose of this paper is to investigate 
the relevant structure of this local solution. 

The overall plan of the paper is as follows. Section 2.1 describes the 
linear inviscid flow produced by the steady upstream distortion field and the 
initial breakdown of the relevant linear 'rapid distortion' theory solution is 
discussed in section 2.2. The appropriate nonlinear, but inviscid, solution 
that eliminates the breakdown is described in section 3. In section 4 we show 
that this solution develops a singularity somewhat further downstream and that 
this singularity can, in turn, be eliminated by a new local solution to the 
problem. 

The viscous boundary layer problem is formulated in section 5. The solu- 
tion to this problem develops a singularity at the position of the inviscid 
singularity and the terminal form of the boundary layer solution is worked out 
in section 5.2. This local similarity solution has a double layer structure 
with an outer inviscid region and a viscous wall layer. A new expansion which 
continues the boundary layer solution through the singularity is worked out in 
sections 5.3 and 5.4. This solution has the same streamwise length scale as 
the external solution and it is shown that the two solutions can be matched in 
an appropriate overlap domain. 

The numerical solutions to the boundary layer problem are described in 
section 6. They show that it exhibits a rapid thickening in the vicinity of 
the singularity, which might be characterized as a kind of "bursting" of the 
boundary layer - a phenomena which has, up to now, been found only in unsteady 
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flows. This might lead to a local transition of the boundary layer and there- 
fore provide a possible mechanism for by-pass transition. These and other 
issues are discussed in section 7. 

2. FORMULATION AND BREAKDOWN OF LINEAR SOLUTION 

We are concerned with the flow over a semi-infinite flat plate of finite 
thickness t* and a leading edge ellipse' whose dimensions are also 0(t*) 
(see fig. 1). The upstream flow is assumed to be nominally uniform, except for 
a small 0(e) steady perturbation, say eU^u^z), in the streamwise velocity 
that depends only on the spanwise coordinate Xz and has a characteristic 
length scale X. u ao (z) is, of course, an order-one quantity, ljfe suppose that 
all lengths have been normalized by X and that the velocity u = {u,v,w} has 
been normalized by the uniform upstream mean flow velocity U^. The flow is 
assumed to be incompressible, with density p, and the pressure p is 
normalized by pU^, The plate thickness t* is taken to be 0(X), and the 
x-coordinate is assumed to be in the streamwise direction with the origin at 
the leading edge while the origin of the y-coordinate is at the flat surface of 
the plate far downstream in the flow. 

Finally, we require that the Reynolds number R\ = U^X/v, where v is the 

kinematic viscosity, be large enough to insure that the viscous effects are 
confined to a narrow boundary layer at the surface the plate that is predomi- 
nantly two dimensional near the forward stagnation point. This will occur if 

In R. << - << R 9 , (2.1) 

which we now assume to be the case. 


2.1. The Linear Solution 

The entire flow is then two dimensional in the vicinity of the leading 
edge, i.e., in the region where x = 0(1), with the three-dimensional effects 
being an 0(e) perturbation of the two-dimensional base flow, say {Uo(x,y), 
Vo(x,y),0}. The viscous effects are confined to a narrow region whose thick- 
ness is o(r; 1/2 Y The solution outside this region should therefore expand 
like V X ' 

U = {U 0 ,V 0 ,0} + e{u 0 ,v 0 ,w 0 ) + e 2 {u 1 ,v 1 ,w 1 } + . . . - {Uq.Vq.O} 

+ eu Q + e 2 u 1 + . . . (2.2) 

and 

p = P Q + ep Q + e Pj + . . . . (2.3) 
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The complex conjugate mean flow velocity £ = IL - iV„ is an analytic 
function of Z - x + i[y +• (t*/X)] that can be expressed in terms of a complex 
potential, say 


W = $ + i* , (2.4) 

where $ is the velocity potential and f is the stream function, in the 
usual way by 


C 


dW 
dZ • 


(2.5) 


For definiteness, we suppose that ? = 0 on the surface of the (body) 
plate and along the stagnation streamline and that $ ■* 0 at the forward stag- 
nation point . 

The first order perturbations are governed by the linearized Euler equa- 
tions and the results obtained by Goldstein (1978) (also see Goldstein, 1979) 
can easily be specialized to the present case to show that the relevant solu- 
tion can be written as 


= v«j) + u oo (z)V A(x,y) , 

(2.6) 

■ -( u o h * v o h> • 

(2.7) 


where 

V 2 (j> = -u ao (z)V 2 A , (2.8) 


A = 4> + 


r* 

1 


1 t 

—00 

U 2 (¥,4») + V 2 (¥,<t>) 


d<t> 


(2.9) 


is t£e Lighthill (1956) - Darwin (1954) drift function and the normal component 
of Ug must vanish at the surface of the plate. 


2.2. Breakdown of the Linear Solution 

The expansion (2.1) is nonuniform in the vicinity of the plate. In fact, 
it is easy to show from these equations, along with the results of Appendix A, 
that the cross flow velocity wq becomes infinite there like 


u m ( z ) 


0 


In ¥ 


as 


¥ ♦ 0 


( 2 . 10 ) 


where a is a constant related to the local potential flow behavior in the 
vicinity of the forward stagnation point and the prime denotes differentiation 
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with respect to z (also see Lighthi 1 1 , 1956). The remaining velocity compo- 
nents, as well as the pressure, remain finite. The strongest nonuniformity 
occurs in the stagnation point region where the mean flow goes to zero. In 
fact, it is easily shown from the results of Appendix A that cpg becomes of 
the same order as Pq at an 0 (e^^) distance from this point. However, the 
linearized equations still give the correct solution (to the order of interest 
here) right up to the edge of the viscous boundary layer and this nonuniformity 
need not be considered further. 

However, another, much more important, nonuniformity develops far down- 
stream in the flow where x becomes large. To understand its structure we 
note that 

* * y* U 0 ♦ 1, * 4 x, A ♦ A + (y) + x as x + « , (2.11) 

where 


- — - - lid* ♦ 4 (y) - I In y as y ♦ 0 (2.12) 

♦ vl / 0 a 


and Aq remains bounded. 

It is now convenient to introduce the new dependent variable 

♦ = ♦ + u oo (z)A + (y) (2.13) 

Equations (2.6) and (2.8) then show that 

V 2 $ = u”A , (2-14) 

T 00 -f- 

and the normal component of u Q will vanish at the plate if 



|^-0 at y - 0 . (2.15) 

9y 

It now follows from (2.12), (2.13), and (2.15) that $ becomes independent of 
x and, consequently, that (see (2.7)) 


-► 



-► 



8<ft 3<fr 
9y’ 3z 


u «\) - 


Pn 4 0 


(2.16) 

(2.17) 


and 
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u ®' 2 

♦ * |3(z) - 2 " In y as y ♦ 0 


(2.18) 


Substituting (2.1) and (2.2) into Euler's equations and using the second member 
of (2.11) along with (2.14) shows that 


3u 


— — + Vd — u ■ Vu — <b ^ — + fd) — u' A ^ ^ — (u , d> , $ — u A 'j 
3x P 1 0 0 ry 3y vz ® +/ 3zJ\ ®’ Y y’ T z ® +/ 


as x + 03 


(2.19) 


It therefore follows from (2.12) and (2.18) that 

,2 


i; < u i + p t ) ♦ - -r ,n y ■ 


(2.20) 


3v 


, 3p, (u„' 2 -uX") , 2 

5 Y ' ny ' 

3. 


( 2 . 21 ) 


3w 3p uV,' 2 

■t “ + •» - r— In y as x * ® and y * 0 

3x 3z 2 


( 2 . 22 ) 


Then since Uj satisfies the continuity equation 


7 • u 1 = 0 , 


(2.23) 


it follows from these results that 


V p. ■* -2 In y as y ■* 0 and x -*■ ® 

1 \ 3 


(2.24) 


and, consequently, that 


♦ d(z) - 


V 'y In y> 


as y -*■ 0 and x -*■ ® . (2.25) 


It now follows from (2.1), (2.2), (2.10), (2.11), and (2.20) to (2.22) that 


2 ,2 

e xu 

u -> 1 + eu (z) - — „ In y (2.26) 
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eu' e xu m u' ’ 9 

w -> — — In y - 5 In y 


a 


a 


(tu'Jy In y > 
p - e 2 d(z) - 


(2.27) 


(2.28) 


as x ■* « , y ■* 0. 

It is clear that this linear expansion must break down when 

-ex In y = 0(1) (2 .29) 

and that a new solution must then be found for this region. 


3. BLENDING LAYER SOLUTION 

To obtain this solution we anticipate that viscous effect will still be 
unimportant and introduce the new gauge functions 

<*(e) = - ln * 0 a S e ■* 0 , (3.1) 

to be specified more precisely below, along with the new scaled variables 

(3.2) 


- _ ex 
x = — 
o 


(3.3) 


n i -o In y - . 

The necessity of using logarithmic variables in problems of this type was 
pointed out by Lagerstrom and Casten (1972) and by Bush (1971) (see also 
Agvawal and Messiter, 1984). 

We simplify the algebra by using the new independent variables 

5 = xn (3.4) 

and n in place of x and n- Then the Euler and continuity equation becomes 

fe a n/a/a La ^ a 1-* fc_ a_ ri/o/a_ A L1_A aJL 

L7 nu af " ave fe n3(i + W 3zJ u " 'La n 3? ’ Van n 35/’ 3zJ P ’ 


and 


3u n^/a l a . aw 

n 35 " ae U + n afr + n = 0 


(3.5) 


(3.6) 
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Equations (2.26) to (2.28) suggest that the solution in this region will be of 
the form 


u = 1 + eu 0 (£,z) , (3.7) 

v = ^e -n/CT nv 0 U,z) , (3.8) 

W = o” n ^0^’ Z ^ ’ ( 3 - 9 ) 


and 


p = e 2 d(z) + | (^-e n/ °) p Q (S,z) . (3.10) 

where ug, vg, wq, and pg are, of course, assumed to remain 0(1) in the 
'blending layer' limit e * 0 with £,n,z held fixed. 

Substituting these into (3.5) and (3.6) and retaining only lowest order 
terms, we obta.n 


and 


9Uj. 

W 


+ w, 


0 3z 


0 



-2 

- 

as 

+ v o 

+ 0 9z 


9w 0 

- S 


as + 

W 0 9z ' 


9 w r 

V_ + a — 

0 3z 


(3.11) 

(3.12) 


(3.13) 


(3.14) 


respectively. Equation (3.13) can be solved for wg, the result can then be 
used in (3.14) to determine vg, which can, in turn, be used in_(3.12) to 
determine pg. Equation (3.11) can, of course, be solved for ug. 


It is well known that the general solution to the 'water wave type' equa- 
tion (3.13) can be written implicitly as 


w 0 - F(z - w q £) 


(3.15) 


where F is an, as yet, undetermined (arbitrary) function of the indicated 
argument. This solution must match onto the linear solution (2.27) in the 
limit as n + 0 with x fixed and the limit x -» 0 with n fixed - both of 
which correspond to the limit £ -> 0. 
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(3.16) 


Expanding (3.15) for small £, we obtain 

w Q = F(z) - F'(z)F(z) n x + . . . 

It therefore follows from (3.2) and (3.3) that (3.9) becomes identical to 
(3.27) when, as yet, the undetermined function F is taken to be 

u ' (z) 

F(z) - - — ^ . (3.17) 

Equation (3.3) shows that n * 1 in the 'small' sublayer 

y = | = 0(1) (3.18) 

of the much thicker region n =* 0(1). Equation (3.15) then becomes 

Wq = F(z - xw Q ) for y = 0(1) . (3.19) 

4. REMOVAL OF SINGULARITIES 

It is well known that the ’water wave' solution (3.15) does not remain 
valid for all £ > 0. In fact, differentiating (3.15), with respect to z, 
shows that 


fjv 

8z 


F' 


1 + £F’ 


(4.1) 


which becomes singular at £ = -1/F' . Here, the prime denotes differentiation 
with respect to the entire argument and we can suppose without loss of gener- 
ality that the initial singularity occurs on the z = 0 plane. This singular- 
ity occurs at the downstream location, say £ s , given by 

£ - Min. g(C) , (4.2) 

S C 


where 


«i = -fin’ ( 

which means that g' should equal zero there (see fig. 2). 

Then, holding £ constant, expanding z in a Taylor series about the 
singular point, and using (4.1) yields 
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where w g is the limiting value of w Q at this point. Then since 


d z _ dz_ 3_ dz g_l n 

j-2 " ,- 3z " F’ u 
dw 0 dw 0 dw 0 


(4.5) 


this, together with (3.14), shows that 


1 „ -2/3 

v 0 * - 3 V 


(4.6) 


and 


*0 * "s + °s z 


1/3 


(4.7) 


as z ■» 0 with £ - 5 S , where |3 S is, of course, an order-one constant. 
On the other hand, (3.13), (3.14), and (4.1) show that 


(4.8) 


and 


w 0 * w s " K s - 5 ’ 


(4.9) 


as £ -> with z held constant. 

We now construct an appropriate local solution that removes this singular- 
ity, or at least delays its formation. The relevant scaling depends on whether 
or not w s is equal to zero. For def ini teness^ we consider only the case, 
corresponding to our numerical example, where w s = 0. The first significantly 
different scaling occurs when 


l 


5 - *s 


0 ( 1 ) . 


(4.10) 


Then equations (4.7) and (4.9), or equivalently (4.6) and (4.8), show that the 
appropriate 'z-scaling' is 


~ Z --Jn = 0(1) ’ 

o 


(4.11) 


and in order to match with (4.7) to (4.9), via (3.7) to (3.10), the dependent 
variables must scale like 


u = 1 + 0(e) 


(4.12) 
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^2 e n/a nv I (C,n.z) 


w = -f 72 nwj( 5 .n,z) 


p = e 2 d(z) + | e n/a j Pj (c ,n , z ) 


where_ vj , wj , and pj remain of order one in the limit as e •» 0 with C. n. 
and z held fixed. Substituting into (3.5) and (3.6) and retaining only the 
lowest order terms now yields 


«, ) a»i _2 - 3 *i 

n V H 11 3z 


P I + 2 n 


*s 9p l 


Z s 8w i _ a *[ 

1 - — v, — - + w. — 1-0 , 

n V 8? 1 di 


£ 3v t 3w r 

- r - + - - 0 - (*- 

1 n 3C 3z 

which clearly generalizes (3.12) to (3.14). The new terms are, in the main, 
due to singularity line curvature effects. 

Equations (4.17) and (4.18) can now be solved for v| and wj and the 
results used in (4.16) to determine pj . Equation (4.18) will automatically 

be satisfied if we introduce the generalized stream function ? by 


f _ Is 3i 

n 3C 


(4.20) 
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Substituting these into (4.17) shows that ¥ is determined by 

(‘ * S* \) * ( f - i* f ?) 4 ( f - H) ( f - H) ' 0 ' <4 ' 21 » 

and in order to match with (4.6) to (4.9) we must require that 

¥ + - r 3 z 1/3 as z * ±<= with £ fixed, (4.22) 

^ s 

and 


¥ + — as £ ■* -<= with z fixed, 


(4.23) 


which clearly represent appropriate asymptotic solutions to (4.21). 


5. THE BOUNDARY LAYER EXPANSION 

Viscous effects must obviously come into play when y becomes suffi- 
ciently small. The boundary layer will initially be two dimensional with the 
cross flow effects producing only a linear perturbation (as in Toomre, 1960) 
until its thickness becomes of the order of the length scale XS introduced 
in the previous section. We therefore suppose that y = 0(1), introduce the 
long streamwise length scale 


and set 


where 


i * _ _ — * 

L " e In 6 



R E 


U L* 


00 

V 


(5.1) 


(5.2) 


(5.3) 


is the 'global' Reynolds number based on L*. It follows that 

6 5 j72 ' (5-4) 

(V >" 1) 

which, in view of (2.1), is consistent with our assumption (3.1) that S << 1. 
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We first consider the region where O^x - £ s ) is greater than 0(a). 


5.1. The Early Three-Dimensional Region 

In this part of the flow we expect the streamwise velocity 

u = U(x,y ,z) (5-5) 

to be 0(1) and in order to satisfy continuity we put 

v = -e S( In S)V (5-6) 

w = -c ( In 5 )W , (5*7) 

where V and W are, of course, assumed to remain 0(1) as e -> 0. Equa- 
tions (3.3), (3.10), and (3.18) suggest that 

p = c^d(z) + 0(e 6 In 5)^ (5.8) 


in this region. 

Substituting these scalings into the Navier-Stokes equations yields the 
three-dimensional zero pressure gradient boundary layer equations 


uiy + v au 

3x 3y 


+ 




„ 3W I. 3W w 3W 

U—+V— +Wg^ 
3x 3y 


2 

31 

3y 2 


(5.9) 

(5.10) 


and 


3U 3V 3W 
— + — + wz 

3x 3y 


(5.11) 


U must clearly go to unity as y ■* 00 in order to match with (3.7). The 
boundary condition for W is a bit more subtle. The known properties of the 
boundary layer solutions suggest that the cross stream derivatives should 
become small as y •+ °°. Equation (5.10) therefore becomes 


3W tt 3W 

— + W 

ax 8z 


(5.12) 


whose solution is given by 


W = G(z - xW) . (5.13) 

This will clearly agree with (3.19) and (5.7) will therefore match with 
(3.9) if we take 
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G = F , 


(5.14/ 


where F is given by (3.17). It is worth noting that the continuity equa- 
tion (5.11) automatically insures that (5.6) will match with (3.8) and (3.14). 

The appropriate boundary conditions for (5.9) to (5.11) are therefore 

U ■» 1, W ■» F(z - xf) as y -» 00 (5.15) 

U = V = W = 0 at y = 0 , (5.16) 

and since the cross flow effects become small as x ■* 0, U must go to the 
Blasius solution in this limit, i.e., 

U ♦ Ug(x,y) as x * 0 . (5.17) 

The solution to this problem should, of course, become singular at the inviscid 
singular point £ s . 


5.2. Terminal Form of Early Three-Dimensional Solution 

It is first necessary to find the terminal asymptotic form of the solu- 
tion to the initial boundary layer problems (5.9) to (5.11) in order to extend 
the boundary layer solution through the singularity. For simplicity, we again 
consider only the case, corresponding to our numerical example, where the cross 
flow velocity vanishes at the singularity. 

We expect the solution to be of similarity form in the vicinity of the 
singular point £ s and the numerical results suggest that it develops a double 
layer structure with a thick, predominantly inviscid, outer region and a some- 
what thinner viscous dominated sublayer. We therefore begin with the outer 
region, put 


y - 



and 


z = 



with y > 0 and seek a solution of the form 


U » f(y,z) , 


(5.18) 


(5.19) 


(5.20) 


V = 


g(y.z) 



(5.21) 
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and 


1/2 


It follows from (5.9) to (5.11) that 


-5 fh + 


3y 


3y 


-yy — + ~z + ? 

9y By 


h(y.z) ■ 

(5.22) 

(h ♦ f If) * - 0 . 

(5.23) 

V 2 ' 3z 


| if) M - 0 , 

(5.24) 


(5.25) 

9z 3z 



The choice of y will, in general, depend on the particular problem, but for 
the case, corresponding to our numerical example, where the singularity lies on 
a symmetry plane, h will behave like h_(y)z and z-derivatives of^ f and g 
will vanish on this plane. It follows from (5.24) that either 9f/3y must 
vanish or 


g = yyf • (5.26) 

Our numerical results suggest the latter, and (5.25) therefore implies that 

h Q + yf = 0 , (5.27) 

which will only be compatible with (5.23) if 

y - 1 , (5.28) 


which, in turn, shows that 


tig = -f on z = 0 . (5.29) 

Equations (4.8) and (4.9) and the inviscid boundary conditions at the plate 
require that 


f + 1 , 


(5.30) 


h -► k(z) , 


(5.31) 


g -► -yh’ (z) 


as y + » and 


(5.32) 
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g~y 2 f ~ y as y ♦ 0 . 
To find the viscous wall layer solution we put 


• («. - 


(5.33) 


(5.34) 


and seek a similarity solution of the form 

, * l+2a 

u - («, - *) F „< 5 ' z > • 
v - Os - • 

, X (3/2)+2a _ „ 

W = (h - V H W< Y ’ Z > • 

It therefore follows from (5.9) to (5.11) and (5.19) to (5.21) that 


(5.35) 

(5.36) 

(5.37) 


2 ~ 9F w / 3 ~ \ 9F W 

-(1 + 2a)F 2 ♦ (G w - aYFJ ~ ^(« w - | *F W ) ^ 


3F 

i w 

3 2 F i 

l 

3z 

“ ~2 
3Y Z 

3H 

3 2 H 

W 

w 

^ * 

s ~2 

9z 

SY* 2 


w 


(5.38) 


-( 2 « * f) F w H » * < C » - “ ?F „> + ( H w * I * F 

and 

3F 9G „ 3F 3H 
-(1 + 2a)F -aY — + -^ + §^^ + -^ = 0 
w 3Y 3Y 3z 3z 

Matching with (5.20) to (5.22) and (5.33) requires that 

F » "if' G » " f2 ■ 

as Y ■* ®, and 

a = 0 . 

The viscous wall conditions require that 

F = G =■ H = 0 at Y =» 0 . 

WWW 


(5.40) 


(5.41) 


(5.42) 


(5.43) 
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It is worth noting that (5.38) to (5.40) have the exact symmetry plane 
solution 

8H 

F - - — = aY ; G = aY^ for z = 0 , 

w 8z 

where a is a constant and that the wall shear must behave like 





x) 


as x -*• £ 

s 


(5.44) 


for all values of z. 


5.3. The Inviscid Inner Solution 

Close to the singularity, where £ - 5s - 0(a), the boundary layer solu- 
tion must match onto the inner solutions (4.12) to (4.15) (rather than onto 
(3.19)) and onto the terminal forms (5.20) to (5.22) and/or (5.35) to (5.37) of 
the boundary layer solutions (5.5) to (5.7). Our preliminary investigation of 
the equations revealed that singularity line curvature continues to play an 
important role even in this relatively thin region. We therefore anticipate 
that the relevant solution will have a double layer structure and that in the 
upper layer will be of the form 


u = Uj(|, yj, z) , 


(5.45) 


v 


2 V j(C. Yj . z) 


(5.46) 


and 


where 


w = 1/2 W i^’ y r 

a 


(5.47) 


p = e 2 d(z) + O^eSj) 2 , 



(5.48) 


and, in order to match with (5.20) to (5.22) we must take 



(5.49) 
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Substituting these into the Euler equations (viscous effects are, of 
course, negligible on the short streamwise length scale (4.4)) we obtain 


and 


U 


. .,(4 

V3y 


I 


5 s 3 \. w 9U I 
- - K + w i — 
yj as; az 


= o , 


u 


a *‘ * »/-! 

as Vay 


i .p 


i y i 


a \ aw, 

J - k + w , ^ 

y T as J 1 1 9z 


= 0 , 


. (j- _ " i.o 

as vay, y, as/ 1 az 


(5.50) 

(5.51) 


(5.52) 


Uj must clearly go to unity in order to match with (4.12). Matching with 
(4.13) and (4.14) can also be enforced if, as we might anticipate from the 
properties of the thin layer equations, Wj becomes independent of yj as 
yj ■* a >. Then Vj will behave like yiV[ j o(CiZ) where Vjq and W[ then 
satisfy 


3W, 3W, 

(1 - V I 0 5 s> -r- * W I - ■ 0 

■ a? az 


and 


(5.53) 


av. 


i o 

r _ P i-ti* + 

1,0 ^S 


as 


aw, 


8z 


= 0 


(5.54) 


as yj ■» <*>, which is clearly compatible with (4.17) and (4.18). The inviscid 
wall condition, of course, requires that 


V, =0 


at 


y r = o , 


(5.55) 


and we anticipate that Uj and Wj will also go to zero there. 

Matching with the upstream solutions (5.20) to (5.22) requires that 


etc . , as S + -® ■ 


U I * 






(5.56) 


(5.57) 
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5.4. The Viscous Wall Layer 


The expansions (5.45) to (5.47) does not, of course, remain valid in the 
region near the wall where y[ •* 0. However, the scalings (4.10) and (4.11) 
should still hold and the upstream solutions (5.35) to (5.37) suggest that u 
and w will be 0(a) and 0(o 1/2 e) there. The viscous terms will then be of 
the same order as the convection terms when y = 0(1) and it follows from the 
requirements of continuity that the wall layer solution should scale like 


u = aUj(C,y,z) , (5.58) 


v = Vj(£,y,z) , (5.59) 

w = eo 1/2 ij(C,y,5) , (5.60) 


Then the Navier-Stokes equations show that the order-one quantities Uj, V j , Wj 
are determined by 


n a_ t v/2- - ^ f- - - -T 

. 1 3? ‘lay y dv 9z V3y y a 5/ J 


{U I ,W I } = 0 , (5.61) 


au, / a L a \ aw. 

a? Vay y a V 3z 


o , 


subject to the viscous wall conditions 


(5.62) 


Uj = Vj - Wj = 0 


at y = 0 , 


(5.63) 


and matching with (5.45) to (5.47) requires that Uj and Wj behave like 


Wj.Uj ~ y as 


y ■* oo , (5.64) 

Matching with the upstream solutions (5.35) to (5.37) requires that, 

(5.65) 


U, - -IK 


*( y ’ (_C) 3/2 ) ’ 
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(5.66) 


Wj - (- 1 ) 


3/2 hA 

W V K) 


3/2 


etc. , as 

Figure 3 summarizes the overall asymptotic structure of the problem devel- 
oped in this and the preceding sections. 


6. NUMERICAL SOLUTION OF THE BOUNDARY LAYER PROBLEM 


The three-dimensional boundary layer problems (5.9) to (5.11) with (5.15) 
to (5.17) was solved numerically using the Keller box method (Keller and 
Cebeci, 1972, Cebeci and Smith, 1974, and Cebeci, Khattab, and Stewartson, 
1981). 

To formulate the numerical problem we first introduce the Blasius variable 

n = y/V^ and write the boundary layer equations as a system of first-order 
differential equations 


P + ep = 


+ »jj! 

, aii 8z ; 


U' = p, 


q’ + eq = 


x^U 


9W i 

0 - 3z 

3x j 


r 


( 6 . 1 ) 

( 6 . 2 ) 


U -/„ 3U w 3W’ 

■J + VS tls > 


(6.3) 


where 1 = 3/3q. 

The boundary conditions are given by equations (5.15) to (5.17) and 

e(x,n = 0 ,z) = 0 . (6.4) 

Note that the definition of e here differs from that used by Cebeci et al . 
(1981) to avoid an apparent singularity in (6.3) as x -> 0. 

A Fourier spectral decomposition (truncated to a finite number of terms) 
is now introduced in the spanwise direction as 

N 

U(x,n,z) - ^ u n (x,n)e inz , (6.5) 

n=-N 

and similarly for the other variables where u_ n = u*, etc. with * indicating 
the complex conjugate. 
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The differential equations for the Fourier components then become 



(6.6) 


(6.7) 


(6.8) 

(6.9) 


The standard Keller box scheme is then applied to the equations for the 
Fourier components in the same way as described by Cebeci et al. (1981) to 
advance the solution downstream from an upstream station where the solution is 
known. To begin the calculation starting conditions are derived by expanding 
the boundary layer solution about the Blasius solution as x •» 0. These 
expansions yield 


U ~ f ' + | nf ' F(z)x 


( 6 . 10 ) 


V ~ - — (f - nf ’ ) + \ (n 2 f ' ' - f - 3nf’)F'(z)\^ (6.11) 

2 V* 


and 


W - f'F(z) 


( 6 . 12 ) 


where f = f ( n ) is the Blasius function. The starting conditions for p, q, 
and e follow from these and (6.2) and (6.3). 

A nonuniform mesh is used in the cross stream coordinate such that the J 
mesh points are distributed according to 

~ ~ vi _ i 

n Q = 0; Hj = hj^ K — , j = l,J . (6.13) 
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This has the effect of increasing the local mesh size from its smallest value 
near the wall, hj, by a fixed percentage (prescribed by the constant K) for 
each subsequent box. 

The introduction of the finite difference approximations yields a 
5(2N + 1) by 5(2N + 1) nonlinear system of equations for the unknown Fourier 
coefficients at each cross stream grid point. The nonlinear system is linear- 
ized about a previous iteration (or initial guess) by Newton's method. The 
difference between two iterations is solved for in an efficient way by the 
block tridiagonal procedure until convergence is achieved. Products are eval- 
uated in physical space and then transformed back to Fourier space using a 
standard Fast Fourier Transform algorithm. 

An advantage to using the Fourier spectral decomposition rather than a 
marching procedure in the spanwise direction is that the need to modify the 
procedure when the cross flow velocity component changes sign (as will happen 
in the cases of interest here) is avoided. The spectral method maximizes com- 
munication in the cross flow variable which then takes on an elliptic nature 
as suggested by Dywer (1981). 


7. RESULTS AND DISCUSSION 

Since vorticity is carried by the fluid particles in an inviscid flow, 
vortex lines that are perpendicular to the plate at upstream infinity become 
infinitely elongated as they pass over the plate-leading to infinite velocities 
and infinitesimal length scales that can only be eliminated by viscous effects. 
This vortex stretching produces the logarithmic singularity in the initial lin- 
ear inviscid solution which, incidentally, is well known in classical rapid 
distortion theory literature (Hunt and Carruthers, 1990). The nonlinear singu- 
larity is new but is also related to this vortex stretching mechanism. It is 
therefore not surprising that its formation is associated with the development 
of small spanwise length scales. 

The solution to the three-dimensional boundary layer problem (5.9) to 
(5.11) and (5.15) to (5.17) depends only on the single function 


which characterizes the scaled upstream disturbance field with the constant a 
being determined by the local potential flow in the vicinity of the forward 
stagnation point. It can easily be scaled out of the problem by replacing e 
by ae in (5.1), (5.4), (5.6), and (5.7). 

Our interest is primarily in periodic-type disturbance fields, and for 
definiteness we choose 


u = cos z 
00 

The initial singularity will then lie at 


(7.2) 
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(7.3) 


= 1 . 


The choice (7.2) has special symmetry properties which may tend to make our 
results less than universal. In particular, it forces the singularity to lie 
on a symmetry plane of the flow which may simplify the terminal asymptotic 
structure of the boundary layer solution. However, these issues are beyond the 
scope of the present work and will be deferred to a subsequent paper. 

The characteristic streamwise length scale L* of the boundary layer flow 
is given by 


L* = - — , (7.4) 

ae In y^R^ea 

when this quantity is renormalized to absorb a into the scaling. The bound- 
ary layer singularity therefore develops on this scale which means that it 
occurs further downstream (relative to the characteristic length scale X of 
the disturbance) as the characteristic amplitude e of the upstream disturb- 
ance and the leading edge/disturbance scale Reynolds number R\ decrease. R\ 
must, of course, always remain large in order for our analysis to be valid. 

The effect of the leading edge geometry is accounted for by the factor a. 

The upstream boundary layer flow is predominantly two dimensional and 
there is an overlap region (in which x << 1 and x >> 0) where it is given by 
the Blasius solution. This solution is therefore the appropriate upstream 
boundary condition for the x-scale boundary layer problem under discussion. 

The computed streamwise velocity profiles are shown in figure 4. They clearly 
show the order-one change in the mean profile shape produced by the small 
cross flow velocity in the external stream. Benney (1984) was among the first 
to show that asymptotically small modulations occurring on sufficiently short 
length scales can lead to order-one changes in the mean flow. However, 

Benney 's modulations arose from interacting Tollmein-Schlichting waves and are 
fundamentally different from the steady upstream distortion effects being con- 
sidered herein. 

The velocity profiles of figure 4 also point to the anticipated localized 
thickening of the boundary layer in the vicinity of the symmetry plane as 
x -*• 1 and to the associated reduction in spanwise length scale. Figure 5 is a 
plot of the upwash velocity profiles. The dramatic local increase in this vel- 
ocity component might be described as a kind of bursting phenomena' which, up 
to now, has only been found in unsteady flows. This might lead to a local 
transition of the boundary layer and therefore provide a possible mechanism for 
by-pass transition. The answer to this question depends on the stability of 
the solution to the local rescaled boundary layer problems (5.41) to (5.48) and 
(5.54) to (5.59). The final resolution of this issue is therefore beyond the 
scope of the present work and will be pursued in a forthcoming paper. 

The cross flow velocity vectors are plotted in figure 6. They show that 
the large upwash velocity is produced by a spanwise flow into the symmetry 
plane. 
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The boundary layer solution develops a double layer structure and the ter- 
minal form of this solution is given by (5.20) to (5.25) in the upper inviscid 
region whose thickness increases like 


1 - x 

and bv (5.35) to (5.40) in the lower viscous region, whose thickness remains 
0 ( 1 ) as 5 * 1 . 

This solution becomes particularly simple on the symmetry plane where it 
implies that 


U = f (y ,0) , (7.6) 

dW _ f(y,0) (7 7 ) 


and 

v _ yf(y.o) ( 7 - 8 > 

(1 - x) 

where the function f is indeterminate and has to be found from the numerical 
solution. The solution in the lower layer is still given by these formulas but 
with f equal to its asymptotic form 

f = ay , (7.9) 

where a is a constant. In this sense the lower layer is no longer a distinct 
asymptotic region but is part of the upper inviscid region to the order of 
approximation of the analysis. This, of course, only holds on the symmetry 
plane . 

To obtain an accurate check on our terminal asymptotic solution we worked 
out the numerical solution to the symmetry plane boundary layer equations 
starting from the upstream Blasius solution. These equations can, of course, 
be solved independently of the flow in the rest of the boundary layer and are 
much simpler than the full three-dimensional equations. It turns out that 
(7.7) and (7.8) constitute an exact solution to these full equations and still 
satisfy the correct free-stream boundary condition, so that the numerical prob- 
lem is especially simple in this case. 

Figure 7(a) is a plot of the symmetry plane streamwise velocity profiles 
(at various values of x) as a function of the composite singularity coordinate 

y(l - x)/Vx. It merely shows that the profiles go smoothly from their initial 
Blasius to their final asymptotic forms. Part (b) of this figure shows that 
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these profiles look quite different when plotted against the unsealed coordi- 
nate y and again reflects the dramatic thickening of the boundary layer. The 
predicted linear behavior at the wall is also evident from this figure. The 
calculated wall shear stress is shown in figure 8. The asymptotic linear 
behavior and the ultimate vanishing of the wall shear, predicted by (5.44), is 
clearly verified by these results. 

The vanishing of the wall shear at x = 1 is indicative of the formation 
of a singularity in the boundary layer solution at this point. But since this 
singularity is directly related to the singularity in the outer inviscid solu- 
tion, it is basically inviscid in nature and it is therefore not clear whether 
or not it is removable. The answer to this question depends on whether the 
local rescaled inviscid problems (4.12) to (4.18) and the corresponding local 
rescaled viscous problems (5.41) to (5.48) and (5.54) to (5.59) possess self 
consistent nonsingular solutions. This constitutes a very difficult numerical 
task that is beyond the scope of the present work. The nonremovability of the 
singularity would require that the upstream flow field predicted by our analy- 
sis be modified due to some form of boundary layer separation that would have 
to be accounted for in the initial formulation of the inviscid flow problem. 
This separation could, of course, be eliminated by terminating the plate before 
the singularity had a chance to form. But even if this is not done the strong 
local instability of the oncoming boundary layer profiles could promote transi- 
tion upstream of singularity and again modify the flow before the singularity 
can occur. We therefore feel that it makes sense to investigate the stability 
of the upstream boundary layer profiles irrespective of whether or not the sin- 
gularity is removable in the asymptotic sense. 

We consider only the destabilization of the classical Tollmein- 
Schlichting waves, since these are likely to be the most unstable modes and it 
is consistent with our asymptotic approach to use the classical high Reynolds 
number asymptotic results given, for example, by Lin (1946) and Reid (1965, 
pp. 281-282). It is appropriate to use the asymptotic formulas associated with 
the lower branch of the neutral stability curve, since we are primarily inter- 
ested in the initial development of the instability and it is not inconsistent 
with our scaling to assume that the Tol lmein-Schl ichting wavelength is short 
compared to X. Then the mean flow can be treated as locally two dimensional 
for purposes of calculating the instability wave growth rates and since the 
most unstable profiles lie along the symmetry plane we consider only this por- 
tion of the flow. 

The asymptotic instability wave growth rates depend only on the Reynolds 
number 



where L* is given by (7.4), and the local wall shear and not on the detailed 
profile shape. We therefore used the wall shear from our symmetry plane bound- 
ary layer calculations in the classical Tol lmein-Schl ichting wave dispersion 
relations to generate the neutral stability curve shown in figure 9. We also 
show the neutral curve for the Blasius boundary layer instability waves at the 
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same streamwise location. Notice that the neutral stability curve moves well 
upstream of the Blasius curve. 

Finally, we note that the present study bears some resemblance to the 
receptivity analyses of Goldstein (1983, 1985} and Goldstein, Leib, and Cowley 
(1987) in that is involves the internalization of free-stream disturbances with 
an attendent streamwise amplification of the perturbed boundary layer flow, but 
is differs from them in that it only involves changes in the mean flow and 
does not depend on the growth of local instability waves to produce the final 
effect . 
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APPENDIX A 


In this appendix we investigate the singular behavior of the drift 
function in^the vicinity of the forward stagnation point where the complex 
potential W behaves like 


W ~ | Z 2 , ( A— 1 ) 

and a is a real constant that is determined by matching with the outer' 
potential flow. It therefore follows from (2.5) that 

c . -aZ - - I^2ai , ( A— 2 ) 

and, consequently, that 

|c| 2 = 2a yV + f 2 . ( A— 3 ) 

Hence 

d<t> 

Kl 2 " 

It follows that 

A 


d<t> 

U 2 + V 2 
0 0 


1_ 

2a 


d<l> 




.2 * 2 

$ +. ¥ 


h 


i 2 + * 2 


* 


( A— 4 ) 


In t 


as 


f + 0 


when 


$ > 0 


( A— 5 ) 
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Figure 2. — Surface of invisdd singularities. 


Figure 3, — Overview of asymptotic structure. 
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(a) x = 0.4. 


(b) x - 0.5. 



(c) x - 0.6. 



(d) x - 0.7. 


Figure 4. — Streamwise velocity profiles as a function of y and z. 
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(a) x - 0.4. 


(b) x * 0.5. 




(C) x * 0.6 


(d) x * 0.7. 


Figure 6. — Velocity vectors in the crossflow plane. 
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Figure 7.— Stream wise vek>city profiles on the symmetry plane as a function of (a) y{1-x)/ V~x~ , and 
(b) y at various values of "x. 
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Figure 9. — Comparison of the lower branch of the neutral stability curve for the Blasius boundary 
layer (dashed curve) with that obtained using the computed wall shear on the symmetry plane. 
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